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Abstract 

These two lectures concern the Color Glass Condensate and the 
Glasma. These are forms of matter which might be studied in high 
energy hadronic collisions. The Color Glass Condensate is high en- 
ergy density gluonic matter. It constitutes the part of a hadron wave- 
function important for high energy processes. The Glasma is matter 
produced from the Color Glass Condensate in the first instants after 
a collision of two high energy hadrons. Both types of matter 
sociated with coherent fields. The Color Glass Condensate is static 
and related to a hadron wavefunction, where the Glasma is transient 
and evolves quickly after a collision. I present the properties of such 
matter, and some aspects of what is known of their properties. 



1 Introduction 



There is a generic similarity between the physical conditions in the early 
stages of heavy ion collisions, and those of the early universe. Such a col- 
lision is illustrated in Fig. 1 . The Color Glass Condensate describes the 
initial quantum mechanical wavefunctions which initiate the collision. This 
is similar to the initial wavefunction in cosmology. At the moment of the col- 
lision there is the initial singularity, analogous to that of the quantum gravity 
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and inflationary stages of cosmology, where quantum fluctuations are impor- 
tant. These fluctuations eventually evolve into the large scale fluctuations in 
the matter distribution produced in such collisions, analogous to the fluctu- 
ations produced in inflationary cosmology which form the seeds for matter 
fluctuations which eventually develop galaxies and clusters of galaxies. Af- 
ter the singularity there is a Glasma phase which has topological excitations 
analogous to those associated with baryon number violation in electroweak 
theory. The Quark Gluon Plasma is analogous to both the electroweak and 
QCD thermal phases of expansion in cosmology, and the deconfinement tran- 
sition is analogous to that which generates masses for electroweak bosons. 
In cosmology there are a variety of phase transitions which occur at various 
times, corresponding to the confinement-deconfinement transition of QCD. 
These various stages in the evolution of matter produced in high energy colli- 
sions are illustrated in Fig. 2. (The labels of this figure are slightly modified 
from the original due to Steffen Bass.) 

There is a new paradigm: Various forms of matter control the high energy 
limit of QCD. These forms of matter have intrinsically interesting properties, 
as well as properties which once understood have extensions to other areas 
of physics such as cosmology. 

Several fundamental scientific questions are addressed by a proper under- 
standing of the properties of these forms of matter. Among them are 

• What is the high energy limit of QCD? 

• What are the possible forms of high energy density matter? 

• How do quarks and gluons originate in strongly interacting particles? 

It is the purpose of these lectures to describe for you the properties of 
these high energy density forms of matter, how they appear in high energy 
physics phenomenology, and what is known experimentally about these forms 
of matter. For reasons of space, these topics are covered only in the most 
generic terms. 

2 The Initial Wavefunction 

In the leftmost part of Fig. 3, I show the various Fock space components 
for a baryon wavefunction. At low energies, the dominant states for physical 
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Figure 1: A schematic picture of the evolution of matter produced in the 
heavy ion collisions. 
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Event Horizon 
Initial Singularity > Quantum Fluctuations 
» 7 x~ - 0.1 fm/c 



■ Initial Nuclei as CGC -> Coherent, High-density Gluons 



Figure 2: The various stages in the evolution of the matter produced in high 
energy hadronic collisions. 

processes have three quarks and a few gluons. In high energy collisions, 
many particles are produced. These ultimately arise from components of the 
hadron wavefunction which have many gluons in them. 

The number of gluons in a hadron wavefunction is usually measured in 
terms of ratio of the gluon energy to the total energy of a hadron in a frame 
where the hadron has very high energy, 

•E ^ gluon j ^ hadron (1) 

Typically, the lowest energy of any gluon is of order Aqcd ~ 200 MeV The 
minimum value of x is therefore x m i n ~ A/ ' Ehadron Therefore, the high energy 
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limit of the description of a hadron is controlled by the low x part of the gluon 
wavefunction. The number of gluons have been measured as a function of 
x, and the results of this measurement are shown in Fig. 4. [1] The values of 
Q 2 correspond to different resolution scales of a virtual photon probe. The 
main point of the figure is that as x decreases, the number of gluons increase. 
The high energy limit of QCD is therefore described by gluonic states with 
large numbers of gluons in them. It is also known from measurements of cross 
section of strongly interacting particles, that their size grows very slowly with 
increasing energy. Therefore the high energy limit is also the limit where the 
gluon density is large. 

By the uncertainty priincple, these components make a gluon wall of 
longitudinal extent of 1/Aqcd as shown in Fig. 3 The density of gluons in 
the transverse plane grows as the energy increases, as is also shown Fig. 5. 



Figure 3: The leftmost figure on this slide illustrates the wavefunction for 
a gluon as a function of energy. The rightmost figure illustrates the experi- 
mentally measured gluon distribution functions as a function of x. 

These gluons must become very tightly packed together. They form a high 
density, highly coherent condensate of gluons, the Color Glass Condensate 
(CGC).[2]-[7] This coherence arises because the high density implies a high 
occupancy of the quantum mechanical gluonic states. In the high occupancy 
limit, the correspondence principle allows the gluons to be described by a 
classical gluon field. Classical fields are coherent. For example, the effects 
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Figure 4: The distribution of gluons as a function of x as seen in the Hera 
deep inelastic experiments. 

of two positive charges add but the effect of a positive and negative charge 
cancel. Because the typical separation of gluons is small, at high enough 
energy as << 1, and the system is weakly coupled. 

Due to the coherence, the Color Glass Condensate is also strongly inter- 
acting. An example of coherence greatly amplifying a very weak interaction 
is given by gravity: The intrinsic strength of the gravitational force is very 
small, but due to the coherent superposition of forces arising from individual 
particles, it becomes a large force. 

There is a typical momentum scale Q sa t, and gluons with momentum less 
than this scale have maximal phase space density, [8]- [10] 

AN 1 (2) 



(PpTCPrxdy as 

As one goes to higher energies, the saturation momentum increases, as more 
gluons are added to the CGC. This is because once the maximum phase 
space density is achieved, then repulsive gluonic interactions forbid greater 
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Figure 5: The density of gluolns in the transverse plane as a function of 
energy. 

occupation. Vacant states of higher momentum must be filled. 

We now have all the information needed to understand the word Color 
Glass Condensate. Color comes from the color of the gluons which make the 
matter. Condensate arises from the high occupation number of the gluons. 
The word glass is because the gluons at small x are produced by gluons at 
larger versus of x. In the reference frame where the hadron is fast moving, the 
high x gluons are also fast moving. The natural time evolution is Lorentz 
dilated. This dilation is passed on down to the scale of the low x gluons, 
which are then evolving on time scales very large compared to their natural 
time scales. This is the property of a glass. Ordinary glass is a liquid on very 
long time scales, but acts like a solid on shorter time scales. 
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Before proceeding further, it is useful to review the properties of some 
kinematic variables. These are commonly used in the description of high 
energy collisions. We define light cone coordinates as 



± -.(t±z)/V2 (3) 



x 

and light cone momenta as 

p ± = {p° ±p z )/V2 (4) 
The scalar product in these coordinates becomes 

x ■ p = xt • Pt — x + p~ — x~p + (5) 
The uncertainty principle for longitudinal light cone coordinates is 

x ± p T > 1 (6) 

It is also useful to define a proper time and a space-time rapidity, useful for 

t > z as 

r = (7) 

and 

The proper time is invariant under Lorentz boosts along the z axis, and the 
rapidity transforms by the shift of a constant. Note that by the uncertainty 
principle, that a particle will have Mt = \jv\ + M 2 ~ 1/r, and that up to 
constants of order one, 

V = \ ln ip + Ip ) ~ ln(p + /M T ) ~ ln(x + /r) ~ r] (9) 

Thus, all the possible rapidities in momentum space and in coordinate space 
are of the same order of magnitude. 

To understand what a sheet of Colored Glass looks like, we can use light 
cone coordinates. Consider F l± ~ E ± B, F + ~ , and F lj . The last two 
terms are of order one. The terms F l± arise from =p derivatives on the vector 
potential. In the frame where typical particle momenta are large, p + is big 
but p~ is small. Therefore we expect that F l+ is big but F l ~ is small. This 
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requires that to a first approximation only the transverse E and B fields are 
large and that 

E±B±z (10) 

These are the non-Abelian generalization of the Lienard-Wiiechart potentials 
of electrodynamics. The essential complication is that the fields carry a color 
index which is random on the sheet of colored glass. The density of such 
fields is determined by saturating the phase space to be of order I /as This 
corresponds to a physical number density of gluons of 

ldN 1 Ql (ID 



tvR 2 dy as 
In Fig. 6, a sheet of colored glass is shown. 




Figure 6: The color electric and magnetic fields on a sheet of colored glass. 

Limiting fragmentation is a well established phenomenon in high energy 
scattering. It is illustrated with the Phobos experimental data shown in 
Fig. 7. [11]- [12] The total multiplicity of produced particles as a function 
of rapidity is plotted in this figure. The data has been shifted so that the 
rapidity of one of the nuclei is at zero. This means rapidities are measured 
relative to that of one of the nuclei. For several units of rapidity from that 



8 




Figure 7: Limiting fragmentation as illustrated in the Phobos experimental 
data on gold-gold collisions at RHIC. 

of the nuclei, the multiplicity distributions are identical. When one gets 
to the central region, corresponding to the smallest x degrees of freedom 
accessible in the center of mass frame for the nuclear wavefunction, the scaling 
disappears. This is at a smaller value of x as the nuclear energy increases. 

The data has a simple interpretation in terms of a renormalization group. 
In the figure, imagine we treat the high x particles as sources. The fields 
produced by the sources are those degrees of freedom in the central region. 
Now as we go from lower to higher energy, there are more degrees of freedom 
associated with the sources, since the produced fields are localized in a few 
units of rapidity. Clearly as we go from one energy to another, we must 
somehow transform what were fields into sources. This happens naturally 
in a renormalization group formalism. One simply integrates out dynamical 
degrees of freedom and converts them to sources as one increases the energy. 
This means that the parameters of the theory which describe the system of 
sources and fields change as one changes the separation scale between what 
we call large x and small x degrees of freedom. 

Such a theory can be mathematically formulated in terms of a path inte- 
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gral. Let us treat the large x sources by a current along the light cone, 



(12) 



where rj is the space time rapidity 



rj = ln(x jx ) 



(13) 



The source is assumed to be independent of x + . This is because of the 
Lorentz time dilation of the sources, and reflects the glassy nature of the 
theory. The theory is defined by the path integral 



In this path integral, the action for the gluon fields in the presence of the 
source p is S[A, p\. The function F[p] is determined by the renormalization 
group equations. The function F is real, and this is a consequence of the 
glassy nature of the system. The parameter A is a separation scale in longi- 
tudinal momenta. Fields with momenta smaller than A are included in the 
fields A, and those with momenta larger than A are in the sources p. 

In leading order, the path integral is approximated by doing a stationary 
phase approximation where the field A is replaced by its classical value for the 
equations of motion. If we do small fluctuations, these corrections are of the 
order asln(p + /A), where p + is a typical momentum scale where we measure 
the classical field. If p + ~ A, then because as << 1, these corrections are 
small, and the classical approximation is valid. If however, p + » A, then the 
classical approximation may break down. In this case, one has to integrate 
out the quantum fluctuations successively, to evolve the theory down to a 
new scale A', where A' is closer to p + . The integrations and change of scale 
for the effective theory define the renormalization group. Remarkably, one 
can get an explicit form for these renormalization group equations, and show 
that the effective action is precisely described by the above path integral. 

The renormalization group equations are for the function 




(14) 



(15) 



are Hamiltonian equations 



dy 




H[d/dp,p]Z { 







(16) 
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Here y = /n(A'/A) For strong to intermediate strength fields, H is second 
order in derivatives. 

It turns out that the Hamiltonian is purely a non-linear second deriva- 
tive theory. There is precisely zero potential. This means that the solu- 
tion for Z corresponds to diffusion. The wavefunction spreads as y — > 
oo. It is this spreading which is the manifestation of the small x problem: 
The saturation momentum grows forever as y — > oo 

Oftentimes, the source function F is taken to be 

F =^2~ i d V d 2 x T p 2 ( V ,x T ) (17) 

This model has many of the features of the full solution of the renormalization 
group equations for F. 

The solution to the CGC renormalization group equations are universal, 
and is a consequence of the diffusive nature of the evolution. This means that 
as y — > oo, F tends to a universal function, implying that the properties of 
the CGC at arbitrarily high energy are independent of what hadron generated 
it. The CGC is a universal form of matter. 

The saturation momentum scale that appears in F, shown explicitly for 
the Gaussian model, generates an infrared cutoff scale. This saturation mo- 
mentum at high energy can be Qs » ^qcd This means that quantities such 
as the total multiplicity of produced particle become computable using weak 
coupling methods. Typically momentum distributions computed in pertur- 
bation theory have infrared divergences at small px- However, these become 
cutoff at the saturation momentum scale, and become perturbatively com- 
putable. For example, the total multiplicity of gluons produced in a collision 
follows from dimensional reasoning: [13] 

1 dN V- (18) 



irR 2 dy as 

The factor of I /as arises because the multiplicity may be determined by 
solving a classical problem, or alternatively, because the phase space density 
of gluons in the initial state has the same factor. We will describe the classical 
problem in the next lecture about the Glasma. 

The x dependence of the saturation momentum itself can be computed 
using the above renormalization group methods. [14] One finds that the sat- 
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uration momentum grows as 



QUv) = Qle Ky (19) 

where y = ln(l/x). This results is to lowest order in weak coupling. Running 
coupling constant corrections change this somewhat, but for intermediate 
values of x, the x, dependence is not greatly modified. 

The Froissart bound on the total cross section also arises very simple. [15]- 
[17] Imagine that we have an impact parameter profile of the matter in a 
hadron. On general grounds, this should fall off exponentially like e~ 2m "" b 
since the lowest mass iso-singlet channel is due to two pion exchange. On 
the other hand, the number of gluons rises like e Ky . Therefore a probe sees a 
matter distribution 



d 2 bdy 

A probe of some fixed resolution scale will see the edge of this matter dis- 
tribution when the distribution is some fixed number. Therefore, the edge 
is determined by b ~ K,y/2m 7T Since the maximal y for beam energy E is 
y ~ ln(E/E ), we see that the cross section saturates the Froissart bound, 

a ~ b 2 ~ ln 2 (E/E ) (21) 

Consider deep inelastic scattering of a virtual photon from a hadronic 
target. The cross section for this should scale like a 7 * p ~ E(Q 2 /Q 2 at ) This 
scaling relation assumes the only scale that determines the physics is the satu- 
ration momentum. There is no added x dependence because the local theory 
which describes the process only has Q sat as a parameter. The geometric 
scaling of Hera data for this process is shown in Dig. 8. [18]- [20] This scaling 
works so lang as x < 1CT 2 . It is also remarkable that the scaling works to 
Q 2 » Q 2 at - The upper limit can be shown to be of order Q 2 ~ Qiatl ^ 2 qcd-> 
and introduces the concept of an extended scaling region. 

The theory of the CGC has provided both a rich phenomenology as well 
as first principles understanding within QCD for the high energy limit. Some 
of the successes of this theory include [7] 

• Scaling properties in electron-hadron scattering. 

• Elastic and almost elastic electron-hadron scattering. 



12 




Figure 8: Geometric scaling in the Hera data.. 

Nuclear size dependence of quark and gluon distributions. 
Distributions of produced particles in hadron and nuclear collisions. 
Scaling properties of hadron-hadron collisions. 

Long range momentum correlations for particles produced in hadron- 
hadron collisions. 

Intuitive understanding of high energy limit for cross sections. 

The origin of the multiparticle excitations which control high energy 
scattering (Pomeron, Reggeon, Odderon). 
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This is a very active area of research. It has as one of its intellectual goals, the 
unification of the description of all strong interaction processes at high energy 
and as such involves different phenomena: electron-hadron, hadron-hadron, 
hadron-nucleus and nucleus-nucleus collisions. 

3 The Initial Singularity and the Glasma 

Before the collision of two hadrons, two sheets of Colored Glass approach one 
another. [21]-[23] Because the phase space density of gluons is large, the gluons 
can be treated as classical fields. The fields are Lorentz boosted Coulomb 
fields, that is Lienard-Wiechart potentials, which are static in the transverse 
plane of the hadrons and have E _L B _L z, where z is the direction of motion. 




Figure 9: After the collision of two sheets of colored glass, longitudinal elec- 
tric and magnetic fields associated with the Glasma form. 

When the classical equations for the evolution of these fields is solved, in a 
very short time scale of order At ~ -^-e~ K l as , the fields change from purely 
transverse to purely longitudinal. This is because the collisions generate color 
electric and magnetic monopole charge densities of opposite sign on the two 
sheets. This description is reminiscent of the flux tube models of color electric 
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fields which are used in phenomenological descriptions of low energy strong 
interactions. At high energies, there are both color electric and magnetic 
field because the electric and magnetic fields were of equal magnitude in the 
CGC, and because of the electric-magnetic duality of QCD. 

To understand how such fields arise, let us first consider the mathematics 
of the one nucleus problem. Assume the sources of the color field are of the 
form 

J+ = 5(x-)p(x T ) (22) 

A solution of the classical field equations is for A ± = 0, and the transverse 
vector potential is a pure two dimensional gauge transform of vacuum field, 
but a different gauge transformation on opposite sides of the sheet at x~ = 0, 

A j = ei-x^-U^xr^UlixT) + 6(x-)-U 2 (x T )V ] Ui(x T ) (23) 

This configuration will generate the color electric and magnetic fields de- 
scribed above. 

In the collision problem, the fields in the backward light cone are deter- 
mined according to Fig. 10. We have chosen to work in a gauge where as 
t — > — oo, the fields vanish. It is easily checked that when crossing the back- 
wards light cone, one gets the proper discontinuity. In the forward light cone, 
one can check that a solution of the equation of motion may be chosen to be 
independent of the space-time rapidity rj. Infinitesmally on the forward light 
cone, the vector potential is A = A 1 + A 2 . This manes that a color electric 
and color magnetic charge is induced on the sheets in the forward light cone 
corresponding to f a b c A{ b E2 C and f a b c A\ b B2 C on sheet 2, and with the labels 1 
and 2 reversed on sheet 1. 

Fields with a non-zero E ■ B carry a topological charge. In QCD, they 
are associated with anomalous mass generation and chiral symmetry viola- 
tion. In electroweak theory, such fields may be responsible for generating the 
baryon asymmetry of the universe. In QCD, they may generate the masses 
of those particles which constitute the visible matter of the universe. Each 
field configuration violates CP. An experimental discovery of the effects of 
such fields would be of great importance. Theoretical ideas for experimen- 
tal signatures are sketchy. [30] These initial longitudinal fields at the time of 
production evolve into transverse fields and eventually form a Quark Gluon 
Plasma. The matter from production until the formation of the Quark Gluon 
Plasma is called the Glasma. [21]- [23] 
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Figure 10: The fields which describe high energy hadronic collisions. 

The initial Glasma fields are unstable against forming rapidity dependent 
fluctuation, and after a time scale of order 1/Q S , instabilities may begin to 
become of the order of the original classical fields. [24]- [28] The origin of these 
Weibel instabilities was originally found for plasmas close to thermal equi- 
librium by Mrowczinski. Quantum fluctuations in the original wavefunction 
can grow by these instabilities, and eventually overwhelm the longitudinal 
electric and magnetic Glasma fields. Perhaps these fields form a chaotic or 
turbulent liquid which might thermalize and isotropize the system. [29] 

The amplification of quantum fluctuations to macroscopic magnitude is 
reminiscent of inflation in the early universe. These quantum fluctuations 
expand to size scale larger than the event horizon during inflation, and are 
imprinted into the fabric of space-time. At much later times, the event 
horizon size scale becomes of the order of galactic size scales. Ultimately, 
these fluctuations drive gravitationally unstable modes which form galaxies 
and clusters of galaxies. 
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These rapidity dependent fluctuations, although they may become large, 
are in the language of cosmology, sub-horizon, and are modified by the late 
expansion. In addition to these rapidity dependent quantum fluctuations , 
there may be rapidity independent fluctuations which are by definition super- 
horizon. These are associated with the initial state of the collision and have 
their origin in the long range Glasma fields. These would be largely unaf- 
fected by subsequent expansion. These may be responsible for the observed 
forward-backward correlations and perhaps the ridge seen in the two particle 
correlation. [31]-[32] 

In heavy ion collisions, analogous fluctuations in the hadronic wavefunc- 
tion might appear as momentum dependent fluctuations. Such fluctuations 
might be frozen into the final state distribution of particles. 
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